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Abstract.  Infrared  Space  Observatory  (ISO)  Long  Wavelength 
Spectrometer  (LWS)  and  ground-based  Fabry-Perot  imaging 
spectroscopic  observations  are  used  to  penetrate  the  extinction 
to  the  powerful  burst  of  star  formation  that  has  occurred  in  the 
extranuclear  molecular  cloud  complex  in  the  galaxy  overlap 
region  of  the  galaxies  NGC  4038/39  (“The  Antennae").  Param¬ 
eters  of  the  starburst  and  typical  molecular  cloud  core  charac¬ 
teristics  are  derived.  It  is  found  that  the  starburst  can  power  the 
infrared  luminosity  of  this  galaxy  system. 

Key  words:  infrared:  galaxies  -  galaxies:  interactions  -  galax¬ 
ies:  starburst  -  galaxies:  evolution  -  galaxies:  nuclei  -  galaxies: 
individual  NGC  4038/39 

1.  Introduction 

One  of  the  most  important  discoveries  resulting  from  IRAS  was 
that  many  galaxies  emit  a  large  fraction  of  their  energy  in  the 
far  infrared  and  that  a  significant  number  of  these  have  lumi¬ 
nosities  of  1012  L0  or  more  (Soifer  et  al.  1987).  Many  of  these 
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galaxies  appear  to  be  colliding  or  merging  pairs.  Theories  on 
the  role  played  by  mergers  in  galaxy  evolution  are  consistent 
with  a  variety  of  origins  for  the  high  infrared  luminosity  ob¬ 
served  in  mergers,  e.g.  shocks  produced  by  cloud-cloud  colli¬ 
sions,  exceptional  bursts  of  star  formation,  or  reradiation  by  dust 
surrounding  a  hidden  active  galactic  nucleus.  In  all  of  these  sce¬ 
narios,  the  presence  of  copious  amounts  of  gas  and  dust  renders 
attempts  to  determine  the  dominant  physical  processes  and  the 
embedded  stellar  content  a  challenging  endeavor.  Far  infrared 
spectroscopy  is  a  tool  well  suited  to  this  task. 

In  this  paper  we  present  ISO/LWS  observations  of  NGC 
4038/39,  a  pair  of  galaxies  thought  to  be  in  the  early  stages 
of  merging.  We  also  present  ground-based  near  infrared  Fabry- 
Perot  imaging  spectroscopy  of  this  system  in  order  to  help  de¬ 
termine  the  spatial  origin  and  excitation  of  the  lines  detected  by 
the  LWS. 

2.  The  Observations  and  Results 

The  ISO/LWS  (Clegg  et  al.  1996)  was  operated  in  its  low  res¬ 
olution  (R  =  A/AA  ~  200)  mode  to  obtain  spectral  line  scans 
of  seven  far  infrared  forbidden  fine  structure  transitions  on  19 
January  1996.  ISO  was  pointed  at  the  position  R.A.(1950)=1  lh 
59m  20(0  Dec(1950)=-18°35'51"  so  that  the  LWS  beam  (Swin¬ 
yard  et  al.  1996)  encompassed  the  extranuclear  molecular  cloud 
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Fig.  1.  Spectra  of  emission  lines  detected  by  LWS. 


complex  in  the  galaxy  overlap  region  (Stanford  et  al.  1990)  as 
well  as  the  nuclei  of  both  galaxies.  Because  data  are  simulta¬ 
neously  obtained  for  all  ten  detector  bands,  almost  the  entire 
LWS  range  was  covered.  The  calibration  accuracy  of  individual 
detectors  was  estimated  to  be  ±30%  based  on  overlapping  data 
in  neighboring  detectors.  An  attempt  was  made  to  correct  for 
these  systematic  errors  by  fitting  the  line  free  parts  of  the  contin¬ 
uum  with  a  second-order  polynomial  in  the  log  F>  -  log  A  plane. 
The  sub-spectra  of  individual  detectors  were  then  scaled  so  that 
their  averages  lay  on  this  polynomial.  Five  of  the  seven  lines 
observed  were  detected.  The  line  spectra  are  shown  in  Figure 
1 .  The  derived  line  fluxes,  their  statistical  uncertainties  (small 
compared  with  the  calibration  accuracy),  and  3cr  upper  limits 
for  undetected  lines  are  given  in  Table  1 . 

High  spatial  resolution  Fabry-Perot  2  pm  spectroscopy  was 
obtained  using  the  NRL/NASM  Fabry-Perot  spectrometer  inter¬ 
faced  to  the  NASA  IRTF  camera  NSFCAM.  The  pixel  size  was 
0.3"  per  pixel  and  the  spectral  resolution  was  R  ~  800.  The  Br7 
and  H2  v=l-0  S(l)  line  images  of  the  nuclear  region  of  NGC 
4039  and  the  extranuclear  molecular  cloud  complex  are  shown 
in  Figures  2  and  3  respectively.  Integrated  intensity  contours  of 
CO(l-O)  from  the  aperture  synthesis  observations  of  Stanford  et 
al.  (1990)  are  superposed  on  the  two  line  images  for  reference. 
The  images  were  registered  with  respect  to  the  CO  interferome¬ 
try  using  the  infrared  position  of  the  nucleus:  R. A. (1950)=  1  lh 
59m  19si7  Dec(1950)=-18°36'27"(±2-3"  )  (S.A.  Stanford,  priv. 
comm.).  Clumps  of  Br7  emission  and  a  diffuse  level  of  line 
emission  are  located  in  and  around  the  extranuclear  region.  The 
brightest  of  the  Br7  sources  are  associated  with  each  of  the  four 
velocity  components  that  make  up  the  extranuclear  cloud  and 


Table  1.  Emission  Line  Fluxes 


Emission  line 

Line  Flux 

10-19  W  cm-2 

Statistical  Error  (1(7) 
10“2°  Wcni”2 

[OIII]  52  pm 

4.9 

2.9 

[NIII]  57  pm 

1.6 

0.58 

[OI]  63  pm 

5.2 

1.1 

[OIII]  88  pm 

4.7 

1.3 

[Nil]  122  pm 

<  0.44  (3<t) 

[OI]  146  pm 

<  0.25  (3(7) 

[CII]  158  pm 

3.7 

0.51 

Br7  2.16  pm 

0.67 

H2  1-0  S(l)  2.12  pm 

0.12 

with  Ha  and  radio  knots  (Bushouse  1987;  Hummel  &  van  der 
Hulst  1986).  Faint  clumps  of  Hi  v=l-0  S(l)  line  emission  are 
coincident  with  three  of  the  Br7  sources.  The  Br7  and  S(l)  line 
fluxes  in  a  30"  aperture  including  the  NGC  4039  nucleus  and 
sources  to  the  northeast  are  given  in  Table  1.  Approximately 
70%  of  this  emission  is  diffuse. 

3.  Discussion 

3.1.  The  Ionized  Gas  and  Starburst  Parameters 

The  presence  of  strong  Br7  emission  in  and  around  the  ex¬ 
tranuclear  cloud  provides  evidence  that  a  localized  starburst  has 
occurred  in  this  region.  The  [NIII],  [OIII],  [Nil],  and  Br7  obser¬ 
vations  can  be  used  to  constrain  the  age  of  the  burst,  assuming 
solar  abundances.  The  CLOUDY  computer  code  (Ferland  1993) 
was  run  with  a  single  star  model  to  assign  characteristic  ionized 
gas  properties  and  an  effective  stellar  temperature  for  the  stars 
responsible  for  the  ionization.  The  observed  line  fluxes  are  fit¬ 
ted,  within  ±2.5(7,  by  a  single  star  effective  temperature  in  the 
range  (3.9  -  4.0)  x  104  K.  The  derived  electron  density  and 
temperature  ranges,  ne  =  60  -  400  cm-3  and  Te  =  7000  ±  1000 
K,  imply  an  HII  region  pressure  range  of  PHii/k=  (0.4  -  3)  x  106 
K  cm-3. 

To  derive  the  parameters  of  the  starburst,  an  evolution¬ 
ary  synthesis  code  (Leitherer  &  Heckman  1995)  was  used  to 
compute  the  starburst  spectral  energy  distribution  for  input  to 
CLOUDY.  For  an  instantaneous  burst  with  a  Salpeter  IMF  ( a= - 
2.35)  for  masses  0.1  -  100  Me,  line  fluxes  were  calculated  as  a 
function  of  the  age  of  the  starburst. 

In  Figure  4  the  logarithm  of  the  ratio  of  the  Br7  to  [OIII]  52 
pm  line  is  plotted  as  a  function  of  the  burst  age  to  illustrate  the 
dependence  of  the  derived  burst  age  on  the  observed  line  ratio 
and  extinction.  For  the  compact  Br7  sources  in  the  extranuclear 
cloud,  we  derive  Av  =  1  -  4  from  both  the  Br-'/Ho  ratio  using 
Ha  images  from  Bushouse  et  al.  ( 1987)  and  the  measured  versus 
predicted  Ha  fluxes  from  the  radio  results  of  Hummel  &  van 
der  Hulst  (1986),  while  Kunze  et  al.  (1996)  derive  Av  =  70  for 
a  model  with  mixed  gas  and  dust  from  Bro/Br  L  Figure  4  illus¬ 
trates  that  because  strong  [OIII]  line  emission  is  limited  to  the 
most  massive  and  short  lived  stars,  the  age  of  an  instantaneous 
burst  is  well  constrained  by  the  Bi'7/[OIII]52pm  line  ratio  to  7  - 
8  x  106  yr,  even  allowing  for  a  large  range  in  the  extinction  and 
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Fig.  2.  Br7  line  Fabry-Perot  image  with  CO(l-O)  aperture  synthesis 
contours  from  Stanford  et  al.  (1990). 

some  Br7  contribution  from  outside  of  our  Fabry-Perot  image. 
The  other  LWS  line  ratios  are  also  consistent  with  this  value. 
Although  we  find  that  a  starburst  with  continuous  star  forma¬ 
tion  with  the  same  IMF  is  not  consistent  with  the  observed  line 
fluxes,  a  finite  rather  than  instantaneous  duration  burst  could 
help  explain  the  age  of  the  youngest  clusters  (age  ~  4  x  106  yr) 
indicated  by  the  Br7  equivalent  widths  of  some  of  the  compact 
Br7  sources  (Fischer  et  al.  1996)  and  by  the  blue  colors  seen 
by  Whitmore  &  Schweizer  (1995).  Based  on  the  far  infrared 
line  fluxes  and  an  adopted  distance  of  26  Mpc  (Vlsr=1550  km 
s-1  and  Ho=60  km  s-1  Mpc-1)  the  total  stellar  mass  produced 
by  the  instantaneous  burst  is  (2.5  ±  1.5)  x  108  M©  and  its 
luminosity  is  (7  ±4)  x  1010  L0.  This  is  about  equal  to  the 
far  infrared  luminosity  (to  within  the  uncertainties)  in  the  LWS 
beam  (see  below),  indicating  that  the  far  infrared  luminosity  is 
powered  by  star  formation  in  this  system. 

3.2.  The  Dust  and  Warm,  Neutral  Gas 

The  LWS  continuum  was  obtained  by  averaging  data  from  all 
line  scans  and  the  associated  multiplexed  data.  The  continuum 
is  best  fit  with  a  i/1  emissivity  law  and  a  dust  temperature  T,/=4 1 
K,  a  typical  value  for  starburst  galaxies  (Stacey  et  al.  1991).  The 
derived  total  IR  luminosity  within  the  LWS  beam  is  3.9  x  1010 
Lq  .  The  derived  temperature  and  emissivity  law  agree  well  with 
that  of  Soifer  et  al.  (1987).  The  total  luminosity  is  lower  than 
that  found  by  Soifer  et  al.  due  to  the  smaller  distance  adopted 
here  and  the  smaller  flux  within  the  LWS  aperture  compared 
with  the  larger  total  IRAS  flux. 

Since  the  peak  wavelength  of  the  observed  [CII]  line  is  at 
the  expected  redshifted  value  for  the  Antennae,  we  conclude 


Fig.  3.  H2  1-0  S(l)  line  Fabry-Perot  image  with  CO(l-O)  aperture  syn¬ 
thesis  contours  from  Stanford  et  al.  (1990). 


Fig.  4.  Log  (Br7  line/[OIII]  52  fim  line)  as  a  function  of  starburst 
age  for  the  instantaneous  burst  discussed  in  the  text.  The  dashed  line 
indicates  the  measured  line  ratio  from  Table  1. 

that  the  observed  emission  is  not  significantly  contaminated  by 
galactic  emission.  The  near  unity  value  of  the  [CII]  158  /t.m/[OI] 
63  jim  ratio  and  the  high  [OI]  63  //m/[OI]  146  //m  ratio  (>  19) 
provide  strong  evidence  that  the  [CII]  and  [OI]  lines  arise  in 
photodissociation  regions  (PDRs).  We  derived  the  PDR  physical 
conditions  and  typical  cloud  properties  following  the  analysis 
of  Wolfire  et  al.  (1990)  (see  Table  2),  using  a  molecular  mass 
of  1.84  x  109  Mq  for  the  extranuclear  region  (Stanford  et  al. 
1990,  scaled  to  our  adopted  distance).  The  range  of  pressures 
derived  for  the  PDRs  are  consistent  with  the  pressure  range 
derived  for  the  HII  regions,  indicating  possible  near  pressure 
balance  between  the  HII  regions  and  the  neutral  medium.  The 
derived  typical  parameters  of  the  warm,  neutral  medium  in  the 
Antennae  are  similar  to  other  starburst  galaxies  (eg.  M82,  Lord 
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et  al.  1996).  The  Antennae  system  falls  in  the  range  of  moderate 
starburst  galaxies  on  the  plot  of  Stacey  et  al.  (1991)  with  an 
overall  [CII]/CO  line  intensity  value  of  2500.  The  young  OB 
stars  must  be  in  close  proximity  to  the  clouds  since  from  Table 
2  the  UV  field  incident  on  the  clouds  is  high  compared  with  the 
average  UV  field,  consistent  with  recent  star  formation  activity. 

3.3.  Excitation  of  the  if  Line  Emission 

For  moderately  dense  (n  <  104  cm-3)  clouds  exposed  to  UV 
photons  from  nearby  OB  stars  where  G0  <  104,  the  FF  v=l- 
0  S(l)  line  and  far  infrared  continuum  emission  are  well  cor¬ 
related  (Luhman  &  Jaffe  1996).  If  the  FF  emission  in  NGC 
4038/39  arises  in  UV  exposed  gas  (PDRs),  then  applying  this 
correlation  to  the  ISO  far  infrared  flux  (scaled  to  the  derived 
area  filling  factor  in  Table  2)  yields  [Is(i)]pz3/;  =  a  few  x  10-5 
erg  s-1  cm-2  sr-1.  The  brightest  FF  emission  is  found  at  the 
NGC  4039  nucleus  (first  detected  at  both  nuclei  in  large  aper¬ 
ture  observations  by  Campbell  &  Willner  1989).  The  surface 
brightness  there  is  Is(i)  =  8  x  10-5  erg  s-1  cm-2  sr-1,  four 
times  higher  than  the  PDR  relationship  would  predict.  In  con¬ 
trast,  in  the  Br7  image  only  a  faint  unresolved  source  is  seen  ~ 
3"northwest  of  the  NGC  4039  nucleus.  We  therefore  argue  that 
in  the  nucleus  of  NGC  4039,  the  emission  is  probably  due  to 
shocks.  These  shocks  are  most  likely  C-shocks  since  the  mea¬ 
sured  S(l)  and  Br7  flux  levels  at  the  nucleus  are  inconsistent 
with  J-shock  models  (Hollenbach  &  McKee  1989).  Does  this 
indicate  that  some  of  the  [OI]  63  /im  flux  in  the  LWS  beam  is 
excited  by  shocks?  Probably  not.  Based  on  the  observed  (un¬ 
corrected)  S(l)  line  strength,  C-shock  models  indicate  that  for 
most  shock  velocities,  [OI]  63  /im  emission  would  be  small 
compared  with  the  flux  level  observed  with  the  LWS  (Draine  et 
al.  1983). 

The  faint  peaks  of  S(l)  line  emission  in  and  around  the 
extranuclear  cloud  region  have  surface  brightnesses  I.sy  i  )  ~  2 
-  4  x  10-5  erg  s-1  cm-2  sr-1,  consistent  with  UV  excitation. 
We  argue  therefore  that  these  FL  emission  regions  and  perhaps 
others  more  deeply  embedded  in  the  extranuclear  cloud  complex 
are  the  likely  sites  of  the  [CII]  and  [OI]  PDR  emission  detected 


by  the  LWS.  These  results  are  also  consistent  with  the  location 
of  the  Br7  sources  in  the  extranuclear  region. 
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Table  2.  PDR  and  typical  cloud  properties 


Physical  Condition  or  Size 

Derived  Values 

npcR 

(0.25  -  1.5)  x  104  cm-3 

G0,  UV  field  flux 

200  -  2500  x  local  ISRF 

T. PDR 

200  ±  60  K 

Ppofl/k 

(0.35  -  3.9)  x  106  K  cm-3 

Mpbr 

8.4  x  107  Mq 

>f>A,  area  filling  factor 

0.04  ±  0.01 

<f>y,  volume  filling  factor 

1.8  x  10“5 

tcore,  molecular  core  radius 

2.32  pc 

Thickness  of  PDR  surface  layer 

0.03  pc 

Ncioud,  number  of  clouds 

1.8  x  105 

Mc ,  cloud  core  mass 

104  M0  (per  cloud) 

<G0>  isrf,  average  ISRF 

45  x  local  ISRF 

